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ABSTRACT 

We present the analysis of two Suzaku observations of GX 301-2 at two orbital phases after the periastron 
passage. Variations in the column density of the line-of-sight absorber are observed, consistent with accretion 
from a clumpy wind. In addition to a CRSF, multiple fluorescence emission lines were detected in both ob- 
servations. The variations in the pulse profiles and the CRSF throughout the pulse phase have a signature of a 
magnetic dipole field. Using a simple dipole model we calculated the expected magnetic field values for differ- 
ent pulse phases and were able to extract a set of geometrical angles, loosely constraining the dipole geometry 
in the neutron star. From the variation of the CRSF width and energy, we found a geometrical solution for the 
dipole, making the inclination consistent with previously published values. 

Subject headings: X-rays: stars — X-rays: binaries — stars: pulsars: individual (GX 301-2) — stars: magnetic 
fields 
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1. INTRODUCTION 

The High Mass X-ray Binary (HMXB) system GX 301-2 
was discov ered in 1969 April during a balloon experiment 
(iLewin et a l. 1971; McClintock, Ricker & Lewin 1971). The 
system consists of an accreting neutron star (NS) fed by the 
surrounding ste llar wind of the B type emiss ion line compan- 
ion Wray 977 (iJones. Chetin & Lillerill974 . A recent lumi- 
nosity estimate de rived from atmospheric mo dels puts its dis- 
tance at ^ 3kpc jKaper. van der Meer & Na iarro 2006), the 
value utilized in this p aper. T he orbital period was estab- 
lished to be 41 days (IWhite . Mason & Sanford 1978) us- 
ing Ariel 5 observations and was refined with the Burst And 
Transient Source Experi ment (BATSE) l o 41.5 days with 
an ecce ntricity of 0.46 (iKoh et alJll997h . iDoroshenko et alJ 
(120 lOa ^ discussed a possible orbital evolution and deter- 
mined an orbital period of 41.482 zfc 0.001 d, ass uming no 
change in orbital period. iKaper. van der Meer & Naiarrol 
( I2OO6) determined that the mass of the companion was in the 
range 39Mo <M < 53 Mq and the radius of Wray 977 was 
R* ~62R0, obtained by fitting atmosphere models. 

The X-ray flux is highly variable throughout an individ- 
ual binary orbit but follows a distinct pattern when averaged 
over multiple orbits (see Figure [T]i- Shortly before the peri- 
astron passage, the X-ray luminosity increases drastically in 
the energy band above ^ 5 keV, as seen in Rossi X-r ay Tim- 
ing Explorer (RXTE)/ All Sky Monitor (ASM) data (iLeahyl 
120021) . The NS passes closest to the companion at a distance 
of ~ O.IR* (Eiydo et al. 1995). Shortly after the periasti-on 
passage, 0orb ~.0^ the X-ray luminosity dips for a short pe- 
riod of time . iLe^ivl ( 12002 ) demonstrated that neither a simple 
spherical wind model nor a circumstellar disk model around 
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Orbital Phase 

Fig. 1 .— RXTE/ASM 1.5-12 keV light curve folded with the orbital period 
of 41.48 days. The two lines indicate the times of the orbital phase when 
the Suzaku observations were performed. The profile shows the strong flux 
increase during the pre-periastrion flare at orbital phase ^ 0.9. For clarity, 
two binary orbits are shown. 



Wray 977 are sufficient to describe the observed variations in 
the folded RXTE/ASM data. An additional stream compo- 
nent is able to account for the sudden increase in X-ray lumi- 
nosity, as the NS passes trough the stream shortly before peri- 
astron and accretes more material. This model also explains a 
slightly higher X-ray luminosity around (poih ~ 0.5, when the 
NS passes through the accretion stream a second time. 

Pulsations with a p eriod of ~ 700 s w ere discovered in the 
Ariel-5 observations (IWhite et al.lll976l) . making GX 301-2 
one of the slowest known pulsars. The pulse period has varied 
drastically throughout the last ^ 20 years (iPravdo & GhoshI 
l2001l;lEvangehsta et alJIOlO) . Prior to 1984, the pulse period 
stayed relatively constant at 695 s-700 s and then spun up be- 
tween 1985 and 1990 to - 675 s. From 1993 until the begin- 
ning of 2008, the change in the spin reversed again, showing a 
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decline. Fermi /Galactic Burst Monitor (GBM) dat^have re- 
vealed that GX 301-2 experienced another spin reversal and 
briefly spun-up with the pulse period decreasing from ^ 687 s 
to ~ 681 s between May 2008 and October 2010. Since Oc- 
tober 2010, the pulse period has shown only small variations 
around ^ 681 s. 

Most recentlv. lGogii^. Krevkenbohm & BeUonil (1201 Ih dis- 
covered a peculiar Iks dip in the luminosity of GX 301-2, 
where the pulsations disappeared for one spin cycle during 
the dip. Several such dips have been previously observed 
in Vela X-1 (Ki-eykenbohm et al. 1999, 2008), where it is as- 
sumed that the accretion on the NS was interrupted for a short 
period of time. 

The pulse phase average spectrum of GX 301-2 is de- 
scribed using a power law with a high energy cutoff. The 
continuum does not show a strong variation in the in- 
trinsic parameters (TjEcut, and /Sfoid) throughout the orbit 
(iMukheriee & Pau H 12004), as seen in two data sets from 
RXTE, taken in 1996 and 2000, sampling most phases of the 
binary orbit. One of the major characteristics of the X-ray 
spectrum of GX 301-2 is the high and strongly variable col- 
umn density of its line-of-sight absorber throughout the orbit, 
indicative of a clumpy stellar wind (A^h = 10^^- lO^'* cm"^). In 
addition to the high column density, a very bright Fe Ka emis- 
sion line can be observed. This line has shown a strong corre- 
lation with the observed luminosity, indicating that the line is 

Siroduced by local clumpy matter surrounding the neutron star 
Mukherjee & Paul 2004). Kreykenbohm et al. (2004) used 
the RXTE data set from 2000 to perform phase resolved spec- 
troscopy and showed that an absorbed and partially covered 
pulsar continuum (power law with Fermi-Dirac cutoff) as well 
as a reflected and absorbed pulsar continuum were consistent 
with the data. 

A cyclotron resonance scattering feature (CRSF) at ~ 
35keV was fi r st dis covered with Ginga (Mihara 1 9951) . 
lOrlandini et al.l (120001) found systematic deviations from a 
power law continuum at ^ 20 and ^ 40 ke V in BeppoSAX, 
where the former could not be confirme d as a CRSF due to 
the proximity of the continuum cutoff. iKrevkenbohm et al] 
(^004) excluded the existence of a CRSF at ~ 20keV and 
showed that the CRSF centroid energy varies between 30- 
38keV over the pulse rotation of the NS. Furthermore, they 
showed that the CRSF centroid energy and width are corre- 
lated. 

We report on two observations of GX 301-2 performed 
with the Suzaku satellite in mid 2008 and early 2009. This 
paper is structured as follows: Section 2 discusses the obser- 
vations and data reduction. Section 3 shows the phase aver- 
aged results. Section 4 discusses the pulse profiles and phase 
resolved spectra. Sections 5 and 6 discuss the results and con- 
clusions, respectively. 

2. OBSERVATION AND DATA REDUCTION 

Suzaku observed GX 301-2 on 2008 August 25 with an ex- 
posure time of - lOks (ObsID 403044010; hereafter Obs. 1). 
The observation was cut short by a set of target of opportunity 
observations and was continued on 2009 January 5, acquir- 
ing an additional ^ 60 ks exposure time (ObsID 403044020; 
Obs. 2). Both main instruments, the X-ray Imaging Spec- 
ti-ometer (XIS; Mitsuda et al. 2007) and the Hai'd X-ray De- 
tector (HXD; lTakahashi et al..,2007) were used in these obser- 
vations. The two observations correspond to orbital phases 



of 0.19 (Obs. 1) and 0.38 (Obs. 2), where Obs. 1 fafls 
into the lowest flux part of the binary orbit (see Figure [T]i. 
The 2- lOkeV absorbed flux was 1.6 x 10~"^ergcm"^ s"' for 
Obs. 1 and - 8 X 10""' erg cm'^ s"' for Obs. 2 (see Table 1 for 
details). Both observations were performed using the HXD 
nominal pointing to enhance the sensitivity of the HXD de- 
tectors. 

The XIS detectors consist of two front illuminated (FI) 
CCDs XIS and 3, and one back ifluminated (BI) CCD, XIS 
1. All three instruments are sensitive between ^ 0.5 keV and 
^11 keV, although the two FI cameras have a higher sensi- 
tivity above ^ 2 keV, and the BI chip is more sensitive below 
~ 2 keV. To minimize possible pile-up, the XIS instruments 
were operated with the 1/4 window option with a readout time 
of 2 s. Data were taken in both 3x3 and 5x5 editing modes, 
which were extracted individually with the Suzaku FTOOLS 
version 16 as part of HEASOFT 6.9. The unfiltered XIS data 
were reprocessed with the most recent calibration files avail- 
able and then screened with the standard selection criteria as 
described by the Suzaku ABC guid^H The response matrices 
(RMFs) and effective areas (ARFs) were weighted according 
to the exposure times of the different editing modes. The XIS 
data were then grouped with the number of channels per en- 
ergy bin corresponding to the half width half maximum of the 
spectral resolution, i.e. grouped by 8, 12, 14, 16, 18, 20, and 
22 channels starting at 0.5, 1,2, 3, 4, 5, 6, and 7keV, respec- 
tively (M. A. Nowak, 2010, private communication). The XIS 
spectral data were used in the energy range of 2- lOkeV for 
reasons described below. 

The HXD consists of two non-imaging instruments: the 
PIN silicon diodes (sensitive between 12-60keV) and the 
GSO/BGO phoswich counters (sensitive above ^ 30keV). To 
determine the PIN background, the Suzaku HXD team pro- 
vides the tuned PIN Non X-ray Background (NXB) for each 
individual observation. In addition, the cosmic X-ray Back- 
ground (CXB) was simulated following the example in the 
ABC guide and both backgrounds were added together The 
PIN data were grouped by a factor of 5 below and by a fac- 
tor of 10 above 50keV for Obs. 1. In Obs. 2 the channels 
were grouped by 3 throughout the whole energy range. GSO 
data were extracted and binned following the Suzaku ABC 
guide. The PIN data energy range of 15-60 keV and the GSO 
data energy range of 50-90 keV were used for the phase aver- 
aged analysis. Due to the short exposure times, no GSO data 
were extracted in Obs. 1 and in the phase resolved analysis of 
Obs. 2. 

2.1. XIS Responses 

GX 301-2 is a source with very large and highly variable 
photoelectric absorption in the line of sight as well as a very 
strong Fe Ka emission line. For spectral modeling, the cal- 
culated responses of the XIS detectors showed a 'leakage' 
emission level (Matsumoto et al. 2006), which stemmed from 
the instrument characteristics. That is, a small fraction of an 
event charge cloud is registered at energies below the peak en- 
ergy, creating a Tow energy tail' in the spectrum. This 'tail' 
is about three orders of magnitude weaker than the intensity 
of the peak count rate of the measured photon. The response 
to the Fe line is therefore described as a Gaussian plus a con- 
stant at lower energies, approximately three orders of magni- 
tude below the peak value and extending to below 1 keV. Due 
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Fig. 2. — (a) Best fit models of the Suzaku GX301— 2 second observation 
for the 3 individual XIS instruments, and the two HXD instruments. At lower 
energies, one clearly sees the constant level in the modeled flux discussed in 
the text, (b) shows that the residuals below ~ 2 keV are very pronounced. 

to the strong Fe Ka line and the very strong absorption in 
this source, the constant level response of the Fe line extends 
above the measured flux below 2 keV. This results in signifi- 
cant residuals below ^ 2 keV, which seem to be stronger for 
the Bl spectrum then for that of the Fl CCD. Figure |2] shows 
the phase averaged spectrum from the second observation of 
GX 301-2 and the best fit models using the HEASOFT gen- 
erated XIS responses as detailed in section 3. 

Together with the Suzaku Guest Observer Facility (GOF), 
which provided for us an experimental response parameter file 
where the 'low energy tail' is ignored, we studied this behav- 
ior in more detail. This new response for all three instruments 
was applied to the individual XIS spectra. Due to the fact 
that the new response completely excluded the missing 'tail', 
we now observed an excess of the data with respect to the 
model below ^ 2keV. The introduction of additional power 
law and/or black body components could improve the fits at 
lower energies. Due to the fact that the true level of the tail is 
not known, these additional components cannot be interpreted 
physically, and we decided to use the original response matrix 
and confine the selected XIS energy range to above 2 ke V. 

3. PHASE AVERAGED SPECTRUM 

3.1. Spectral modeling 

Phase averaged broad band spectra in the 2-60 ke V (Obs. 1 ) 
and 2-90 ke V (Obs. 2) energy ranges were obtained. From the 
technical description of the XIS instrument^ it is known that 
small discrepancies, e.g., in fitted power law slope, have been 
observed between Fl and Bl XIS instruments. A difference 
of ^ 0.05 in the power law index has been observed in cali- 
bration data and was also observed here. When modeling the 
Fl and Bl instruments with a common power law index, the 
residuals of the Bl instrument show a systematic deviation. 
However, only the Fl or Bl instruments could be modeled to- 
gether with the HXD instruments. Due to this fact and the 
higher sensitivity of the Fl XIS instruments above 2 keV, we 
concentrated our discussion on the results obtained with the 
two Fl instruments. Each data set, Fl and Bl, was modeled 



individually with the HXD data to compare the differences in 
the continuum. We found that when using the Fl and Bl in- 
struments individually with the HXD instruments, in both ob- 
servations the best fit spectral values are consistent with each 
other within error bars, indicating that the usage of the HXD 
reduced the observed discrepancies in the power law parame- 
ters. 

The continuum model for pulsars can so far only be mod- 
eled with empirical models, consisting of a powerlaw with 
a cutoff at higher energies, which is typical for this kind of 
source (Coburn 2001). Three empirical models are widely 
used: the simple high energ y cutoff (h iqhecut) and the 
Fermi-Dirac cutoff (f dcut, iTanakalll9 86) both are used in 
combination with a simple power law component. In addi- 
tion, the negative-positive exponential powerlaw model NPEX 
dMihara 1995) is a slightly more complicated model including 
the power law component. For the data analyzed in this work, 
the best results have been obtained with the f dcut model: 



IPL 



exp(£'-£'cut/£'foid)+l' 



(1) 



where ApL is the normalization at 1 keV, E^m is the cutoff en- 
ergy, and Efoifj is the folding energy of the Fermi-Dirac cut- 
off. For the Obs. 1 and Obs. 2 Fl data sets the best fit x''s 
were 260 and 542 with 245 and 267 degrees of freedom (dof), 
respectively (see Figure [3] and |4]i. Replacing the smoother 
Fermi-Dirac cutoff with the highecut model in Obs. 2,re- 
sulted in slightly worse residuals (Ax^ w 15) compared to the 
best fit results of the FD-cutoff model (Figure [3] d) and an 
emission line-like residual at ^ 35 keV, which is most likely 
due to the sharp break of the power law at the cutoff energy 
in this model. In comparison, the NPEX model shows even 
worse residuals (Ax^ ~ 200) and only results in reasonable 
fits when the exponential curvature at higher energies is inde- 
pendent in the partially and fully covered component. 

The best fit cutoff energy is very close to that of the ob- 
served CRSF feature (see below) and these two parameters 
are rather strongly correlated (see Figure |5j. To avoid a de- 
generacy of the CRSF values due to a changing cutoff energy, 
Ecut was frozen in all the fits to the best fit value from the 
Obs. 2 Fl spectrum (29.2 keV). 

Previous observations of GX 301-2 indicated the existence 
of clumps in the stella r wind ( Kreykenbohm et al. 2004t 
iMukheriee & PauJ 120041) . which were modeled using par- 
tial covering absorption in addition to fully covered pho- 
toelectric absorption of the smooth stellar wind. In the 
present analysis, the low energy portion of the spectrum 
also required a partial covering component (Figure |3}' and 
1^), which was modeled using the TBnew model (Wilms 
et al., 2011, in pre p.fl an updated version o f the existing 
TBabs model ( Wil msTAUen & McCrayll2000t) . In addition, 
a non-relativistic, optically-thin Compton scattering compo- 
nent cabs was included, as is necessary for column densi- 
ties A^H > 5 X 10^^ cm"^, where the plasma becomes Comp- 
ton thick and part of the emission is scattered out of the line 
of sight. The A^h values from TBabs and from the cabs 
component were set equal and treated as one model compo- 
nent for the smooth stellar wind (A^h.i = TBabs 1* cabs 1) 
and a second component for the clumpy partial coverer (A^h.2 
= TBabs2*cabs2). The A^h,i and A^h,2 column densities 
were left independent of each other. For spectral fitting the 
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Fig. 3. — Data and best fit model for Obs. 1. a) shows the data and the 
best fit model. The inset shows the residual for the different emission lines 
observed in XIS data, where the Fe Ka line is already included in the model, 
b) shows the best fit residuals without partial covering and the CRSF. c) has 
only the CRSF added, and d) has both componenst included. 

abundan ces of|W ilms. Allen & McCrayl (2000) and the cross- 
sections (I Verner et al. 1996) were used in all data sets. 

Both observations show residuals to the continuum in the 
35-40keV energy range, which we interpret as t he pre- 
vious ly observed CRSF (Mihara 1995; Krevkenboh m et al] 
12004"). Modeling these residuals with an absorption line with 
a Gaussian optical depth (gabs) improved the residuals sig- 
nificantly in both observations (see Figures[3]and|4]l. The best 
fit values of £cRSF = 42.5!^;^ keV and 35.2!i;4keV for Obs. 1 
and Obs. 2, respectively, are consistent within 90% confi- 
dence intervals. The widths of the CRSFs, ctcrsf = 9.5i!^2^ 

keV 

(Obs. 1) and 7.8;!;()9keV (Obs. 2), are also consistent within 
errors. An increase of ctcrsf with hig her centroid energy, 
as previously observed in RXTE data dKrevkenbohm et al] 
|2()04), could not be detected, although the best fit values hint 
at such a behavior. Table [T] shows the best fit values for the 
continuum with a gabs CRSF line where the given errors are 
90% confidence values. To calculate the significance of the 
CRSF in the first and weaker observation, the null hypothe- 
sis approach was applied, where 10,000 spectra were created 
with Monte Carlo simulations using the best fit parameters 
without the CRSF. Gaussian uncertainties were used for the 
individual model parameter. Each model was fitted with and 
without the CRSF line to compare how much an inclusion of 
the line actually improves the fit. In ^ 99.6% of all fits, no 
CRSF feature were observed with a larger optical depth than 
the observed lower limits of r = 10.3, making the existence 
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Fig. 4. — Same as Fisure [3]for Obs. 2. In this case b) shows the residuals 
without partial covering and the CRSF, c) has only the CRSF added and d) 
has both components included. 
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Fig. 5. — Contour plots of ficRSF vs Ecm- The different contours indicate 
the 1, 2, and 3 sigma contours for the FI spectrum of Obs. 2. The X indicates 
the best fit values for both parameters: Ecut = 29.2 keV and £cRSF = 35.2 keV. 

of the feature significant with a ^ 3. Including the CRSF in 
the real data, the best fit improved by a A^^ = 100. A similar 
improvement of was only observed in ^ 0.1% (a ^ 3.3) of 
the simulated spectra, concluding that the observed line in the 
first observation is indeed real. For the second observation, 
the line was even more pronounced, where the exclusion of 
the CRSF increased the by 200. 

The CRSF can alternatively be described with the Loren- 
zian shaped eye labs XSPEC model CMihara et al...l990.) . 
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TABLE 1 

Phase averaged spectral parameters with the gabs model component for best fit spectrum 



Parameter 


403044010 (lOks) 


403044020 (60ks) 




FI 


BI 


FI 


BI 




16.6!^ 


16.6 (frozen) 


20.9+3| 


20.9 (frozen) 


Abund Ca 


1.55 (frozen) 


1.55 (frozen) 


1 55+0-39 
• -^-0.18 


1 98+0.41 
'•'-0.49 


Abund Fe 


1.17 (frozen) 


1.17 (frozen) 


1 ,7+0.04 
'•"-0.03 


1 04+0-03 

'•" -0.04 


A'h,2[1022 cm-2] 


76.9+3« 


76.9 (frozen) 


28.4+1-0 


28.4 (frozen) 


r 




n 09+0.09 
^■°^-0.08 


96+0 06 

"•^"-0.04 


98+0 03 
"•'°-0.03 


ApL.l [10^^]' fully covered 




1 OO+0.20 
'-'°-0.20 


1 1+1-1 

'■'-0.5 


4 70+0-91 
^•'"-0.79 


ApL,2 [10^^]' part, covered 


an <:+16.1 
' -10.0 


90.0!^,- 


215.5+25,1 


266.4+17-9 


£cut [keV] 


29.2 (frozen) 


29.2 (frozen) 


29.2+13;5 


29.2 (frozen) 


fifold [keV] 


10 3+3' 


8.3+?'« 




5.6^0-2 


£cRSF [keV] 


42 5+3- 


39.7+3-5 




34.9^;9 


o-CRSF [keV] 


n c+1.6 
^■-'-2.8 


8.6+?'« 




7 4+0.5 
' ■^-0.5 


fCRSF 


on a+17.7 
-19.0 


18.9!'^ 


12.0^V 




£sKa [keV] 


9 o^+O.OI 

-0 02 


9 Ti+o:03 

-0 03 


9 QQ+O.Ol 


2.34t«l 






9 44+0.96 


5-36ti 


4 94^+0.68 


EQWi^^ [eV] 


860!^«* 


1109!^313 


67irf« 


'<JJ-445 


^ArKa [keV] 


'^■'^ -0.03 


-'■"^-0.05 


3 00+^1-01 


3 04+O-02 

-'•"^-0.02 


/ArKc<[10-5]^ 


1 79+0.62 
-0.69 


1 oo+<'-'"* 

-""-0.48 


9 74+0.55 
'^-0.54 


4 <;i+0.94 
••^■^-0.88 


EQWLk, [eV] 


I 7Q+502 
' '^-164 


1 10+3"** 


38+29 


64^38 


EciKa [keV] 


T 71+0.05 
^ - '^-0.04 


3 76+<'-'" 
^■'°-0.06 


3 70+0-02 
-'■'"-0.01 


-1 71+0.03 
■'•"-0.03 


/CaKa[10-5]2 


1 71+0" 
'■"-0.83 


2 45+1-33 


4 79+0.95 
^■'^-0.99 


4 95+1.68 
^•"-1.71 






100+154 
'""-100 


14+10 
'^-8 




EpeKa [keV] 


f. lQQ+0.007 
"■-^'^ -0.006 


"■^'"-0.010 


6 409+0.003 
"■^"^-0.001 


f: 494+O.003 
"■^^^-0.002 


O-peKa [eV] 


13+f 


< 31 


< 5 


< 7 


/FeK«[10-4]2 




6.96^:« 


19 20+0-29 
'^-^"-0.30 


99 9+0.61 
^^■^-0.58 




241+53 




133!27 


142+10 
'^^-10 


£^FcK/3 [keV] 


7 06+" °3 
' -""-0.03 


7 04+<'-'"i 

' -"^-0.04 


7 09+0-'" 

'■"^-0.01 


7 09+O-02 

' -"^-0.02 


/FeK/3[10-*]2 


1 1Q+0.29 
'■''-0.35 




2 65+0-27 


4 50+0-53 
^•-""-0.56 




43+33 


56!i3 


2i;^ 


32+9 


^NiKa [keV] 






7 46+O-02 
' ■^"-0.02 


7 51+0 -11 

'•-"-0.06 


/NiKa[10-^]2 






1 57+0-25 
'■-^ '-0.25 


67+0-62 

"•"'-0.62 


EQW^.K,, [eV] 






15+9 


5 7+13-4 

-'•'-5.7 


Flux^ absorb. 


1 61 +'"'4 
'■°'-0.15 


1 63+<'-03 
'■"-'-0.07 


7.98^-S? 


9 04+0 04 
^•"^-0.06 


Flux* iQj^^y unabsorb. 


16.3+2'0 


i6.o;o-8 


31.8:i;3 


35 9+1-0 


Cxiss/Cpin/Cqjq 


0.93;«;O^/1.18t«l/- 


-/1.12!<J[]^/- 0.951 


+0.003 /1 T1+0.02/1 on+0.12 
-0.003' '- '-0.02/ '••''-0.12 


-/I 17+''-°2/l 27+0-11 
' '-"-0.02/ '-^'-0.11 


xVdof 


260 / 245 


145/130 


542 / 267 


328 / 150 


(l)Units are ph keV"' cm"^ 


s"' , (2) Units are ph cm"~ 


s ' , (3) Values of EQWs are determined relative to the abs. 


continuum. 



(4) Absorbed and unabsorbed flux units are [10 l" erg sec 1 cm 2], (5) Values of C are with respect to XIS for the FI fits and XIS 1 for the BI fits. 



Cyclabs is described by the centroid energy Eqrsf, the 
width (TcRSF, and the resonance depth tcrsf, similar to the 
gabs parameters. The best fit continuum parameters are 
consistent with the best fit values determined with the gabs 
component. The cutoff energy with a best fit value 
of Ecut = 31.2;|;2 9keV is consistent with the value obtained 
with the gabs model. The observed centroid energies of 
£cRSF = 35.0!^:5keV and 31.6!i:°keV are of the order of 
10-20% lower than the energies obtained with the gabs 
model. This discrepancy stems from a different calcula- 
tion of the line centroid energy a nd is described in detail in 
iNakaiima. Mihara & Makishimal (120101) . The width (Tcrsf = 
11.2!];gkeV and 12.3![:;^keV for Obs. 1 and Obs. 2, respec- 
tively, is bigger than with the gabs model. The I dof val- 
ues of 260 / 245 and 548 / 267 for the two observations using 
the cyclabs model and thus were not a better fit when com- 
pared to the gabs model. For the final discussion we use the 
values determined by the gabs model. 

Several emission features were observed in the residuals of 
the fits in both observations and were subsequentially mod- 
eled with Gaussian emission lines (see Figures[3]and|4] inlay). 



The width of each line was set equal to that of the Fe Ka emis- 
sion line, while the intensities were left to vary independently. 
Energies were loosely constrained around the expected values 
of neutral material to avoid runaway of the line energies. 

A constant (const) was applied, taking small instrumental 
differences in the overall flux normalization into account. The 
constant was fixed at 1 for XIS and was left free for the other 
instruments. Cxis3, Cpin and Cqso in Table 1 are the cross 
calibration constants with respect to XIS for FI fits and with 
respect to XIS 1 for the BI fits. 

The final model had the form: const*A^H,i (PLl-i- 
A^H,2*PL2)* f dcut*gabs-i- 6xGaussians, where the 
power law indices of PLl and PL2 are set to be equal to each 
other The individual power law normalizations were inde- 
pendent and were used to calculate the fraction of the partial 
covering. 

Tablem summarizes the best fit values for both observations 
and for the individual FI / BI data sets. For the interpretation 
we will concentrate on the FI data for both observations, for 
reasons mentioned above. 
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3.2. Spectral results 

The best fit spectral parameters were generally consis- 
tent with previous observations, with the exception of the 
cutoff energy. With a value of 29.2;!;2^]^ keV, Ecut was 
significantly high er than the ^ 20keV va lue measured in 
previous RXTE dMukheriee & Paull 12004 and BeppoSAX 
dLa Barbera et al.l 12005) observations, although these obser- 
vations used slightly different spectral models for the cutoff 
energy. .Krevken bohm et al. ( 2004) used also the here ap- 
plied Fermi-Dirac cutoff, resulting in cutoff energies of 10- 
15keV. The folding energy for Obs. 1, EfM = 10.3!i] ikeV, 
is consistent with values obtained for a sim ilar orbital phase 
with BeppoSAX dLa Barbera et all l2005h . In Obs. 2, the 
value of iifoid = S.Ti^^jkeV is consistent with the RXTE 
data for the pre-periastron flare and the periastron passage 
jKrevkenbohm et al.ll2004l) . 

The observed column densities for the absorption due to 
the smooth stellar wind (A^h.i) were consistent between both 
observations. A larger difference could be observed in the 
column density associated with the absorption due to the 
clumped wind (A^H.2), where the best fit value of Obs. 1 is 
almost a factor three larger than that of Obs. 2. In addition 
to the H column density in the TBnew model, the relative 
Ca and Fe abundances were also left independent in Obs. 2. 
The best fit values are measured from the Ca and Fe K edges 
at 4.1 keV and 7.1 keV, respectively, and were slightly higher 
than solar abundances: 1.55^jj [g for Ca and 1.17;!;q q3 for Fe. 
For the first observation, the Ca and Fe abundances could not 
be well constrained and were fixed to 1 .55 for Ca and 1.17 for 
Fe. 

The partial covering fractions can be calculated from the 
measured normalization values of the two power law compo- 
nents, ApLi andApLi: 



Cvr. Frac. 



1PL2 



iPLl 



+ A 



(2) 



PL2 



The factor is almost unity in both observations due to the dom- 
inance of the ApL2 value. Similar large co vering fractions have 
been observed in XMM-Newton data (Ftirst et ani2011 b') and 
in the RXTE data (Mukherjee & Paul 2004), indicating that 
the covering fraction does not change significantly for differ- 
ent parts of the orbit. The spectrum softened slightly from 



r = O.S3-^^"(^ for Obs. 1 to O.ge^U JJ^ for Obs. 2. 

The Fe Ka Hne was detected at 6.409!;j jjjj^ keV (Obs. 2), 
which is, given the instrumental gain systematics, consistent 
with neutral Fe. In addition, the following lines have also 
been observed in both observations: S Ka line, Ar Ka line, 
Ca Ka line, Fe K/3 line, and the Ni Kg Une, all with energie s 
consistent with neutral material (e.g. lKaastra & Mewdll993h . 
Note that the Ni Ka line was not detected in the fainter first 
observation. The observed line intensities are summarized in 
table[T]with 90% confidence errors. The values are consistent 
between the FI and BI instruments. The widths of the lines 
were set to be equal to the Fe Ka width within each observa- 
tion, which was found to be 13^2^ eV for Obs. 1 and had an 
upper limit of 5 eV in Obs. 2 (FI values). The Compton shoul- 
der to the Fe Ka line, as observed by the Chandra observatory 
{Watanabe et al. 2003), was not sig nificantly detec t ed in e ither 
observation. The observation of .Watanabe et al.l (l2003 h was 
performed in the pre-periastron phase, where the luminosity 
was higher than in the two Suzaku observations. A Compton 
shoulder has also been observed with XMM-Newton in data 




Fig. 6. — Pulse profiles for different energy bands for Obs. 2. The two 
peaks are indicated in the top panel. The bottom pannel shows the selection 
of phase bins which were used for phase resolved spectroscopy. 

taken during another pre-periastron flare (iFiirst et alJl201 lb|) . 

4. PHASE RESOLVED ANALYSIS 

Barycentric and binary corrected light curves for differ- 
ent energy bands were extracted for both o bserv ations us- 
ing the orbital parameters from ' Koh et al] (11997 ) with an 
updated orbital pe riod and periastron time to determined by 
iDoroshenko et al.l (12010a). Due to the strong variability of 
the pulse period on short time scales, individually determined 
pulse periods were used for each observation. An accu- 
rate pulse period could be established for Obs. 2, but not 
for Obs. 1 due to its short duration of 10 ks. For the sec- 
ond observation the calculated pulse period had a value of 
P = 685.4 ± 0.9 s. The Fermi/GBM instrument measures the 
pulse period of GX 301-2 on a regular basis, resulting in val- 
ues of P = 687.28 s for Obs. 1 and = 685.75 s for Obs. 2 
(Finger, 2011, priv. comm). To create the pulse profiles, the 
Fermi/GBM-provided pulse periods and epoch times for each 
observation were used. 

We studied the energy dependent pulse profiles for Obs. 1 
and Obs. 2. Both pulse profiles are very similar, showing a 
double peak shape where the main peak (PI) is broader below 
lOkeV. The second peak (P2) stays rather constant in width, 
but increases its relative intensity toward higher energies. As 
an example, the pulse profile for Obs. 2 in different energy 
bands is shown in Figure |6l Comparing these pulse profiles 
with previous RXTE and BeppoSAX data showed that the 
general shape is consistent through all parts of the orbit. 

The data of the longer Obs. 2 were divided into 10 equally 
spaced phase bins and individual spectra were extracted for 
the XIS and PIN instruments. Figure |6] shows the individual 
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phase bins in the lowest panel. This division resulted in an 
exposure time in each phase bin of 5-6ks for each individ- 
ual instrument. No GSO data were used in this analysis due 
to the reduced exposure time per phase bin. For spectral anal- 
ysis, the phase averaged model was applied to all phase bins, 
resulting in best fit values summarized in Tabled as well as 
Figure [T] Again the cutoff energy was frozen to the phase 
averaged value of 29.2 keV to avoid the previously discussed 
degeneracy with /scrsf- 

The two absorbing components, A^h.i from the smooth stel- 
lar wind and A^h,2 from the partial coverer vary between 20 
and 40x 10^^ cm"^, but show an anti-correlated trend through- 
out the first peak. Although the error bars are rather big, one 
can see an indication that A^h.i follows the flux in PI, whereas 
the Nh,2 value dips at the same time. 

In contrast to Kreykenbohm et al. (120041) . the power law in- 
dex r is varying strongly throughout the pulse. The values 
are 0.9-1.3 throughout the first peak and drop suddenly to 
0.6-0.8 for the second peak (Figure|7]l. Measuring the second 
peak to be significantly harder than the first peak supports the 
behavior seen in the pulse profile, where the intensity of P2 in- 
creases at higher energies, becoming similar to the intensity in 
PI. Power law normalizations are relatively small and badly 
constrained for the first power law component. The partially 
covered power law normalization seems to follow the flux, 
showing a higher value throughout the first peak. The folding 
energy does not vary significantly throughout the orbit and 
shows values between 5 and 6 keV (not shown in Figure 7) 

The CRSF energy varies between 3 0-40 keV, similar to 
the v alues observed in the RXTE data (iKrevkenbohm et aLl 
120041) . Table 2 shows that the addition of the CRSF in flie in- 
dividual spectra improved the for most of the phase bins, 
except in phase bin 6, which falls in the gap between both PI 
and P2. Phase bins 1 and 7, also minima in the pulse profile, 
only show a small improvement in the best fit when the CRSF 
was included. The CRSF energy is not very well constrained 
in phase bins 1 and 6 due to the lack of sufficient statistics. 
For the geometrical discussion below, these two phase bins 
are ignored. In all other phase bins the addition of the CRSF 
improves the fit. The energy changes very smoothly through- 
out the pulse and does not follow explicitly the observed flux 
in each phase bin. The CRSF energy increases during the 
main peak and reaches the maximum at the falling flank of 
PI. During P2 the CRSF energy decreases until it reaches a 
minimum at the dip between P2 and PI. 

To estimate the significance of the CRSF in the phase re- 
solved analysis, again the null hypothesis method was applied 
in two out of the 8 remaining phase bins, which are good rep- 
resentatives of all data points. Phase bin 2 is a good example 
for a shallow CRSF, whereas phase bin 5 is an example for a 
phase bin with a higher flux. Similarly as with the phase aver- 
aged data, 10,000 spectra were simulated and modeled with- 
out and with the CRSF. The interpretation of the observed 
feature as being due to stochastic fluctuations with a depth 
similar to the observed values could be ruled out with a prob- 
ability of over 99% {a ^ 3) in both phase bins. This leads to 
the conclusion that the CRSF can be observed with sufficient 
significance to allow the interpretation below. 

5. DISCUSSION 

5.1. Phase averaged continuum 

In both observations very strong absorption can be ob- 
served. A^H.i, which is the absorption due to the smooth stellar 
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Fig. 7. — Phase resolved spectral parameters for Obs. 2. Top panel shows 
the 20 - 40 keV pulse profile, where the individual phase bins are numbered. 
Units of power law normalization are ph keV"' cm"^ s"' . The intensity of 
the Fe line is units of ph cm"^ s"' . For clarity, two pulses are shown. 
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TABLE 2 

Phase resolved spectral parameter from Obs. 2 for FI data. Same model as phase averaged data. 



Parameter PBl PB2 PB3 PB4 PB5 PB6 PB7 PB8 PB9 PBIO 

iVH.i[10-2W] 19.2+j»,2 20.0+^-^ 36.l!:*-5 2(12+11 324+f7^ 22.8+f;^ 253^^ 1Z251 21.4!^-^ 24^6+ff" 

iVH,2[10-2/cm2] 32.6!f!; 30.9!^-;^ 21.3!|3.8 25.0!2i:o l&j+H 26J+l^^ 32Jtii 35.2!^:o 29.5:3| 31.6!2| 

r 1 99+0.08 1 :in+0.b4 n nQl+0-03 1 9n+0.b4 1 07+0-10 n 7A+0"05 n n 67+0 05 1 15+0.06 

i ^-^^-O.U ^•^'^-0.04 "-^^-0.04 "-^^-0.04 ^-^"-0.04 ^■'^'-0.22 "''^-O.OS "■-'°-0.04 "'"'-O.OS ^-^^-0.06 

ApL,i[10-2]'ful.cov. 0.9rlf, 1.35!313 10.12!«« O.irf^^' 8.69!2«^ 1.92r;« 2.03t_\il O-ielJ-?? 0.83!^'?^ 2.36!2.« 

ApL2[10-2]' par. GOV. 24.4+*-^ 44.8^^] 17-4:|-g 27.4!2j2^ 27.7!^-^ 10.8!?-j] 10.6!|-|; 11.3|'/ 19.3!3-0 

£M,[keV] 7.4!i:^i 6.13!!;;i 5.6!«j 5.5:»i 5.9:°;^ 4.1^1 4.5^ S.O^i 4.4^ 5.0^^ 

£cRSF[keV] 38.7+^-3 29.S+[^ 32.1+j'^ 34.6+3-^ 38.o;2j 38.3!f,'' 35.4+*-i 35.7+j-^ 31-i:|o ^'^■^-oi 

acRSF[keV] 10.3+|^ 5.1t[\ 6.4^| 7.3:|t ('■^^l ^^-^^ ^-^-[t ^-^Ti ^J^'i ^-^-ol 

rcRSF 33.4!]5.i 3.5!ii 6.1^^ ^.4^_\l ^3^1 30.5^^ \1 .O^fi 2\.^^,l 11.4!f| 3.9!i| 

/FeK.[10-3]2 1.91^,% 1.82+!;;1!] 1.86ti^ l-S^i? 1.72^-;«^ 1.76!«;|f 1.87:!];i? 2.16^:!! 2.11^;j| 2.07+0:1!] 

FH-.o.eV 4.93!«;M 7.77:!];«« 10.03^:« 11.27!0;3^ 10.93^:3' 7.76:];01 e.56^l\ ^.41^,g 7.62^;>o 4.89^:«5 

Cxisa/CjiN 0.94/1.15 0.94/ 1.15 0.95/1.19 0.94/ 1.21 0.94/1.21 0.94/ 1.19 0.95/ 1.16 0.95/1.14 0.94/1.10 0.94/ 1.11 

x"/dof 233/243 288/243 353/243 312/243 253/243 303/243 283/243 292/243 264/243 247/243 

X^/dof no CRSF 24 1 / 246 307/246 39 1 / 246 341 / 246 273 / 246 305 / 246 29 1 / 246 344 / 246 297 / 246 270 / 246 

1 ^ i '7 i ~\ C\ i 

(1) Units are ph keV" cm , (2) Units are ph cm"~ s" , (3) Units are 10 erg cm sec" , (4) Values of C are with respect to XIS for the Fl fits. 



wind, is consistent within both observations. A^h,2, the absorp- 
tion from the partial coverer, is significantly stronger in the 
first observation, establishing the existence of an additional 
component in the line of sight. Leahy & Kostka (2008) used 
RXTEIkSM and PCA data to show that a stellar wind and a 
stream model component can describe the observed count rate 
throughout the orbit and that the increase of the column den- 
sity at orbital phase ^ 0.2 is mainly due to the stream compo- 
nent. The findings in both Suzaku observations are supportive 
of this theoretical picture. 

GX 301-2 shows v ery strong luminosity va r iability on 
very sh ort time scales dKrevkenbohm et al]|2004t iFiirst et al.l 
1201 Ibl) . where the column density increases to up to 50% in 
periods of lower activity. Both RXTE and XMM-Newton 
data showed a much higher column density in the pre- 
periastron flare than in the two data points observed with 
Suzaku. One possible explanation is that the wind is in- 
deed so variable and clumpy that these variations are just not 
properly predictable, especi ally in times of hig h activity, such 
as the pre-periastron flare. ' Leahy & Kostkal (l2008i) investi- 
gated RXTE/PCA Nh values of the absorbed component from 
archival observations for different parts of the orbit and found 
no increased column density in the pre-periastron flare. Ob- 
servations throughout one full binary orbit could help to un- 
derstand the variations in A^h on time scales of days. 

The observed luminosity dependence of the power law in- 
dex r and the folding energy is similar to other sources, 
such as V0332H-53 (Mowlavi et al. 2006) and 4U0115H-63 
(iTsygankov et al .' 2007). The power law index F shows a 
hardening for the first observation, whereas the fo lding energy 
£foid is slightly higher when compared to Obs. 2. ISoong et al] 
( [1990) observed the Efoid variation in Her X-1 phase resolved 
HEAO-1 data and concluded that the parameter is dependent 
on the viewing angle of the accretion column and can be used 
to describe the plasma temperature of the system. In the case 
of GX 301-2 the smaller Efou value in Obs. 2 could indicate 
a lower plasma temperature which can be interpreted that the 
X-ray emission region is further up in the accretion column 
dBasko & Sunvaevll 19761) (see also Section 5.2). 

The softening of the power law index T with increased lu- 
minosity is also in agreement with the basic model of the ac- 
cretion column that the plasma temperature decreases with in- 
creased height. With increased luminosity, the rate of accreted 



material would increase and the amount of soft photons cre- 
ated by the lateral walls of the relatively taller column would 
increase, leading to a softer spectrum, as is observed. 

For the second observation the abundances for Fe and Ca 
were left independent and a slight overabundance was ob- 
served. _TakingJiito_£onsWera^ that the abundances used 
in Wil ms. Allen & McCravl (l2000l) are derived from the inter- 
stellar medium (ISM) a small overabundance from an evolved 
star may be expected. 

5.2. Variations in the CRSF parameters 

In both observations the CRSF was clearly observed and 
the gabs absorption component improved the overall fit sig- 
nificantly. Strong magnetic fields {B ^ 10'^ G) exist close to 
the NS magnetic poles, where photons at energies close to the 
Landau levels are resonantly scattered from the line of sight 
and result in an absorption line-like feature in the spectrum. 
This feature provides a direct method to measure the magnetic 
field strength close to the NS surface, where the fundamental 
centroid energy can be described as 

£cRSF=11.6keVx^x^ (3) 

where z is the gravitational redshift at the scattering site, with 
z ~ 0.3 for typical NS. Using the ^crsf values obtained in the 
phase averaged spectra, the calculated magnetic field has a 
strength of 4.76!!;:^^ x lO'^G (Obs. 1) and 4.10![!;],^ ^ 10'^ G 
(Obs. 2) which is consistent within errors for both observa- 
tions. 

The CRSF parameters are consistent with previous results 
using RXTE data and are lower than the observed B eppoSAX 
(la B arbera et al.ll2005h and INTEGRAL (Doros henko et aP 
.201 OB) values of 45-50keV. These observations were ob- 
tained during pre-periastron outburst, where the luminos- 
ity of the source is much higher than in the two Suzaku 
observations presented in this paper A luminosity depen- 
dence of the phase-averaged CRSF centroid energy is ob- 
served in multiple other sources, where an anti-correlation be- 
tween CRSF and luminosity w'as observed in V 0332-H5 3 and 
4U 01 15H-63 (ITsygankov et al.ll2006HMihara et alJl2007h. and 
a positi ve correlation wa s obser ved in Her X-1 dStaubert et alj 
120071) . iKlochkov et alj ( 1201 Ih have confirmed, from pulse 
to pulse variability studies, such correlations for V 0332H-53, 
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4U 01 15+63 and Her X-1, and also found an anti-correlation 
for A 0535+26. Although the CRSF values are consistent 
with a single value in both Suzaku observations, a hint of 
an anti-correlation can be seen in the data. Compared with 
the BeppoSAX and INTEGRAL data, the CRSF centroid en- 
ergy might have a posi tive correlation. Paralle l to the CRSF- 
luminosity correlation, iKlochkov et al.l (1201 Ih also observed 
in pulse to pulse data of four sources that the power law index 
r shows an opposite luminosity dependence than the /scrsf- 
In the Suzaku observations discussed here, T is observed to 
soften with increased luminosity, another indication that a 
possible negative CRSF-luminosity correlation exists. 

The two opposite scenarios seem to depend on the luminos- 
ity, i.e., if the observed luminosity is below or above a critical 
luminosity (CL) which can be derived from the Eddington lu- 
minosity of a system dBecker & Wolff 2007). The CL is of 
order ^ 10^^ erg s"', but depends also on the accretion geom- 
etry (spherical or disc) as well as the height and diameter of 
the accretion column. 

Above the CL, the infalling proton density becomes so large 
that the protons begin to interact and decelerate, creating a ra- 
diation pressure dominated shock region above the magnetic 
pole. This region of increased densit y is most likely the region 
where the CRSF is created (Basko & SunyaevI 19731) . With in- 
creasing luminosity, the shock region moves higher up in the 
accretion column, where a smaller local magnetic field value 
results in a lower observed CRSF centroid energy, as is ob- 
served in V 0332+53 or 4U 0115+63. The observed CRSF- 
luminosity dependence as well as the F-luminosity depen- 
dence in GX 301-2 is very similar to the scenario described 
here. 

Below the CL, the accreting matter slows down via hy- 
drodynamical shock, 'Coulomb friction' or nucl ear collision 
(iBasko & SunvaevI [19731: iBraun & Yahell E^S?). With in- 
creasing accretion rate, the luminosity increases and the de- 
celeration region is pushed closer to the NS surface, where 
the magnetic fields are higher. This would result in a posi- 
tive CRSF-luminosity correlation, as observed in Her X-1 and 
A 0535+26. 

With a distance of 3 kpc, the intrinsic unabsorbed 2 - 10 keV 
luminosity of ~ 2 x 10^^ erg s"' is significantly below the typ- 
ical CL of ~ 10^^ erg s"' . Although it has been observed that 
sources with similar luminosit ies can show the opposite cor- 
relation ( K lochkov et al.l201 Ih due to individual critical lumi- 
nosity values, the large luminosity difference compared to CL 
puts this source in the same regime as Her X-1 and A 0535+26 
and at odds with the notion of a well-defined CL. The ob- 
served variation in F and possibly in iicRSF would then stem 
from another mechanism. 

5.3. Emission lines 

The existence of multiple fluorescence emission lines, es- 
pecially at lower energies, where the absorption is very domi- 
nant, indicates that the source of the emission originates from 
a region where the column density is not very large, i.e. the 
outer l ayers of the stellar wind in our line of sight ( Fiir st et al] 
1201 Ibl) . If the emission lines would be embedded deeper in 
the stellar wind, the lines, especially at lower energies would 
have to be significantly stronger, to be detected at all. On the 
other hand, if the line emitting region is in the outer layers 
of the wind, the incident soft X-ray flux would be drastically 
reduced, making the equivalent width very large, as observed 
in the low energy emission lines (Table 1). A possible expla- 



nation would be that the emission region is very large for the 
lines, and maybe spread over the entire surface of the stellar 
wind. 

The most dominant line emission stems from the Fe Ka 
transition at 6.4 keV. The observed energies in both observa- 
tions are consistent with the emission in neutral material. The 
ratio of the intensities of Fe K/3/Fe Ka for both observations, 
0. 16!g:jj^ for Obs. 1 and 0. 14!jj{][ for Obs. 2, is consistent with 
the fluorescing materia l being neutral or only slightly ionized 
(iKaastra & Mewell993l) . The equivalent widths (EQW) of the 
two Fe Ka lines show that the intensity relative to the contin- 
uum is reduced by a factor of ~ 2 in the se c ond o bservation. 
Based on emission line widths, lEndo et al] (120021) estimated 
in ASCA data that the Fe emission o riginates w i thin 10' " cm 
of the continuum emission source. iFurst et al.1 (1201 Ibl) use 
XMM-Newton pulse by pulse observations to observe a cor- 
relation with increased flux of the source, showing that the 
Fe emission region is not far from the X-ray source. In the 
Suzaku observation, however, the Fe Ka line flux does not 
directly follow the observed luminosity. The increase in ab- 
sorbed 2- lOkeV luminosity is a factor of ^ 5, whereas the 
increase in the fluorescence line intensities is only 2-3. This 
difference in intensity change indicates that the distance be- 
tween continuum and fluorescence emission lines is large, 
especially when c ompared to the proposed ^ 0.3 It-s from 
lEndo et aTl ( 120021) . In the 6 months' time between the two 
observations, however, the overall emission geometry could 
have changed, accounting for the difference in the fluores- 
cence intensities. If the emission originates from a confined 
region in the stellar wind, observations taken during two dif- 
ferent orbital phases would each yield a different distance to 
this region. 

In addition to the strong Fe Ka line, many other emission 
lines are observed in the XIS spectra in both observations. 
The observed intensities of these lines show the same behav- 
ior, where the flux for the second observation is ^ 50% or 
less compared to the first observation, although the observed 
uncertainties for the different lines, especially for S Ka and 
Ar Ka, are very large. Furthermore, the S Ka line falls into an 
energy regime where calibrati on uncertain ties have been pre- 
viously observed (e.g.. ISuchy et al.ll201 II) . so that a detailed 
analysis is not possible at this time. 

5.4. Pulse profiles 

Historically, GX 301-2 shows very strong pulse to pulse 
variations and an intensive flaring behavior on very short time 
scales (Kreykenbohmet al. 2004; Fiirst et al. 2011b), e.g., 
throughout the pre-periastron flare. The observations obtained 
here do not show such flaring behavior, and only show the 
regular pulsation throughout the observation. A pulse pe- 
riod could only be established for Obs. 2 of the Suzaku data, 
which was consistent with Fermi/GBM data. With the GBM 
pulse period for the first observation, pulse profiles for both 
observations could be produced for different energy bands. 
The two peaked pulse profiles show a similar behavior for 
both observations, where the second peak is getting strong at 
higher energies. This behavior is i n contrast to many sim - 
ilar sources, such as 4U 0115+63 (Tsvgankov et al.' '20071^ 
V 0332 +53 (iTsygankovet al. 2006). 4 U1909+07 (Furst et d] 
1201 lal) . and 1 Al l 18-61 (ISuchy et al.li201 ih . where the pulse 
profile turns into a single peak profile at higher energies. A 
reduction in the intensity of the second pulse has also been 
observed in BeppoSAX data for different parts of the orbit 
tLa Barbera et al..2005,) . 
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Fig. 8. — Calculated phase resolved magnetic fields from the CRSF energy 
and the best fit for a simple magnetic dipole model. 

5.5. Geometrical constraints using a simple dipole model 

The phase resolved analysis of the second observa- 
tion showed strong variations in several spectral parame- 
ters throughout the pulse profile. As previously observed 
by iKrevkenbohm et al. (2004), the CRSF centroid energy 
varies by ^ 30% where the highest energy is detected at 
the peak and the falling flank of PI. Such behavior has 
been observed in mul tiple other sourc es, e.g.. Vela X-1 
(ILa Barbera et al.i'20()3|K 4U 0352-H309 (ICoburnll200l |). and 
Her X-1 ([Klochkov etd]|2008l) . A model for the CRSF vari- 
ations is based upon the change in the viewing angle through- 
out the pulse and thus different heights of the accretion col- 
umn are probed, yielding a different local magnetic field ob- 
served for each phase bin. 

A simple approach to derive a possible geometry for the 
neutron star and the magnetic field uses the variation in the 
observed magnetic field throughout the pulse phase. In this 
case the very smooth and sinusoidal variation (Figure 8) can 
be modeled as a simple dipole where the total magnetic mo- 
ment /i is calculated from the phase averaged CRSF energy 
value. 

The variation of the CRSF energy is then fitted by chang- 
ing the geometrical angles of the system until a best fit con- 
verges. Appendix A discusses the model in detail, introducing 
the three free parameters: O, the viewing angle between line 
of sight and the NS spin axis, the inclination angle a of the 
magnetic moment with respect to the spin axis, and the an- 
gle /3, indicating the 'lag' of the B-field plane with respect 
to the ephemeris, i.e., the observed shift in pulse phase (see 
also Figure 9) The best fit values for /3 were either -4.6° or 
175.4°, depending on the starting values of the fit (see Ap- 
pendix A), where the latter value is equivalent to the former 
when the dipole polarity is flipped, a and 9 angles showed a 
very strong interdependence, and the best fit 0/a pairs were: 
-14.6°/67°, 14.6°/113°, 67°/- 14.6°, and 113°/14.6°. AH 
best fit values have a w 4 with 7 dof. These pairs do show a 
degeneracy of 0± 14.6° and 90±23°, which can be explained 
with a geometrical symmetry, when calculating the magnetic 
field for each individual phase bin. The first and last two pairs 
can be each treated as the same geometry, rotated by 180 de- 
grees. Figure |8] shows the calculated magnetic field values 
from the CRSF centroid energies and the theoretical values 
for B from the best fit angles for a simple dipole model with 
e = 67°, a = -14.6° and /3 = -4.6°. 

The phase resolved energy and width of the CRSF show 
a strong correlation, where the width varies in phase with 



the CRSF energy (Figure. 7). IKrevkenbohm et all (120041) 
observed a similar correlation, where the magnetic pole 
was observed under different viewing angles $, where $ is 
the angle between the lin e of sight and the magnetic axis. 
iMeszar os & Nagell (Il985h showed that the anisotropic veloc- 
ity field of the electrons in the accretion column leads to a 
fractional Une width of: 

J81n2-^|cos$|=/t|cos$|. (4) 

i'cRSF V »IeC^ 

The (Tcrsf/£^crsf ratios vary in the 0.16 - 0.24 range, cor- 
responding to variations of ±20% throughout the pulse phase. 
The two outliers with large errors in phase bins 1 and 6 were 
ignored. The variation of $ can be estimated as - a and 
Q + a throughout one pulse phase, where the magnetic pole 
rotates around the NS spin axis, tilted by 9 with respect to 
the line of sight. 

Using the average ratio of 0.2, we calculated the values 
of k for the two geometries with 9 = -14.6° (k = 0.21) and 
9 = 67° (k = 0.51). For each geometry, we used the variation 
of the (Tcrsf/^crsf ratio (±0.04) to calculate the variation 
in resulting in values of $ = 15° ± 25° and $ = 67° ± 5°. 
These variations are smaller then the values expected from 
the the geometrical discussion, but show a similar behav- 
ior, where the variation of the angle is larger for the smaller 
value of $. From Equation 4, we were also able to esti- 
mate a possible plasma temperature for the two geometries 
from the calculated values of k. The plasma temperature of 
kT ^ 4keV for 9 = 15° is very similar to the observed fold- 
ing energy Em^, which is an indication of the plasma temper- 
ature (Burderi et al. 2000, and references therein). The value 
of kT ^ 24keV for the 9 = 67° geometry is much larger and 
is not consistent with the observed fifoid- Assuming that the 
NS spin axis is aligned with the inclination of the binary sys- 
tem, a value of 9 = 15° would put the inchnation at ! '--^ 75°, a 
value which is o nly marginally above the upper l imit of 72°, 
as determined bv lKaper. van der Meer & Naiarro (2006). 

The strong Fe line was detected in all 10 phase bins. The 
measured line flux did not change significantly throughout 
the pulse phase, while the 2- lOkeV flux varied by more 
than a factor of 2. From this we conclude that the dis- 
tance to the Fe fluorescence region is greater than ^ 700 It-s 
(^ 2 X 10 '-^ cm), from the N S. This is in conflict with the con- 
clusion of lEndo et alj (|2002[) that the emission region must be 
closer than 10^° cm from the NS, based upon the width of the 
emission lines measured by ASCA. Their assumption is that 
the emission region is close to the NS and that the material 
free falls onto the NS surface, encountering much faster ve- 
locities, which caused the observed broadening of the lines. 
The smaller measured Fe Ka line width in the phase averaged 
spectra is more consistent with the assumption that the line 
broadening is due to the terminal velocity of 300 -400 km of 
the line driven wind (Parkes et al. 1980). 

6. SUMMARY 

We presented results of two Suzaku observations of 
GX 301-2 taken shortly after the periastron passage. The 
spectra were modeled with a partially and fully covered power 
law with a Fermi-Dirac cutoff and a CRSF at ^ 35keV. 
The column density of the smooth stellar component only 
marginally changed and the clumpy wind component signifi- 
cantly changed between both observations. Flux dependen- 
cies of the CRSF and F were not significant, but do hint 
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at similar correlations such as those observed in V 0332+53 
and 4U 01 15+63, although the observed luminosity is signifi- 
cantly below the calculated critical luminosity of this system. 
They are in the realm where one might expect the energy of 
the CRSF to decrease with decreasing flux. The variations in 
the pulse profiles and the CRSF throughout the pulse phase 
have a signature of a dipole magnetic field. Using a simple 
dipole model we calculated the expected magnetic field val- 
ues for different pulse phases and were able to extract a set of 
geometrical angles, loosely constraining the dipole geometry 
in the NS. Model constraints derived using the observed ratio 
of CTcRSF to /icRSF, together with calculated plasma tempera- 



tures, favor a solution wherein the spin axis is tilted ■ 
our line of sight. 
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APPENDIX 
A. GEOMETRICAL DIPOLE MODEL 

The variability of the CRSF energy throughout the pulse phase can be used to understand the geometry of the neutron star. We 
assume a simple dipole magnetic field around a spherical neutron star, however we want to emphasize that these assumptions are 
very rudimentary, and effects such as the displacement of the dipole field from the neutron star center are not taken into account 
here. 

The total magnetic moment /i of the neutron star can be calculated with: 



(Ai)=B-/?^s(VT+3cos2$)= 6.776 X 10^°G cm' 



(Al) 



using the average magnetic field from the phase averaged CRSF centroid energy (B = 4.1 x lO'^G) and a typical neutron star 
radius of /?ns = 10 km. The value of the last part, including the azimuthal angle $ of the magnetic pole was estimated to the 
average value of 1 .54. 

The total average magnetic moment is divided into three components based on the cartesian coordinate system in Figure |9] 



'sin(Q!)cos(/3) 1 
m = ^ I sin(a)sin(/3) 
sin(/3) 



(A2) 



where a is the angle between the dipole field and the spin axis and /3 is the rotation angle from the ephemeris, i.e. the x-axis. 
Next, the vector n was calculated for each individual phase bin used in the phase resolved analysis. 



^sin(e)cos(0)^ 
sin(9)sin((/i) 
sin(</)) 



(A3) 



The individual components of n are the projections of each phase bin on the three axes, as seen from the line of sight of the 
observer. The angle corresponds to the angle between the line of sight and the spin axis. For each of the individual phase bins, 
the magnetic field was calculated using the position vector n indicating the position of the magnetic pole and the vector m, which 
indicates the magnetic moment: 

3 nlm n| m . 

5phase = - —T (A4) 



The calculated magnetic field values were then fitted to the measured magnetic field, varying the three angles O, a, and /? until 
the best fit converged. To avoid wrong solutions due to local minima, the starting parameter of all three angles were varied in 
steps of 10° before the fit was started. The best fit solutions resulted in the values that are discussed in detail in section 5.5. 
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